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We show that the modulation of the phases of the laser beams of ultra-short pulses leads to 
modulation of the two photon fluorescence intensity. The phase modulation technique when used in 
multi-photon microscopy can improve the signal to noise ratio. The technique can also be used in 
multiplexing the signals in the frequency domain in multi-focal raster scanning microscopy. As the 
technique avoids the use of array detectors as well as elaborate spatiotemporal multiplexing schemes 
it provides a convenient means to multi-focal scanning in axial direction. We show examples of such 
uses. Similar methodology can be used in other non-linear scanning microscopies, such as second or 
third harmonic generation microscopy. 


Interaction of matter with coherent light field leads to 
excitation with unique phase signature. Different spec¬ 
troscopic techniques have been developed that rely on the 
interferometry of the wave-packets in the material sys¬ 
tems generated by a sequence of laser pulses with con¬ 
trolled relative phases HHS]. In these experiments the 
systems under investigation are resonantly excited and 
the phases of the laser pulses are manipulated to select 
out certain pathways in the excitations that are of in¬ 
terest laisE]. Similar approach can be taken in non¬ 
resonant excitations, such as multi-photon absorption. 
We use phase modulated laser pulses to induce two pho¬ 
ton absorption in CdSe quantum dots. We show that the 
intensity of the subsequent fluorescence from the quan¬ 
tum dots is modulated at characteristic frequencies. This 
allows us to separate the fluorescence arising from the 
different interaction pathways. We apply the technique 
in multi-photon microscopy (MM) to address the spe¬ 
cific challenges related to simultaneous raster scanning 
of multiple points [Mo]. 

When using phase modulation in MM we split a pulsed 
laser beam from a femtosecond oscillator (Synergy from 
Femtolasers) into two replicas using a 50/50 beam split¬ 
ter (Figure [^. The frequencies of the pulses in the 
two beams are sheared by few tens of MHz by using 
acousto-optic modulators (AOMs). The electric fields 
El {t) and E 2 (t) in a pair of pulses in the two beams can 
be written as Ei{t) = A{t) cos{{uj + and E 2 {t) = 

A{t — t)cos{{uj + — r)), respectively, where A{t) 

is the envelope of the electric field, uj is the carrier fre¬ 
quency, r is the time delay between the pulses, and (pi 
and p 2 are the frequency shifts imparted by the modula¬ 
tion of the phase. The time delay between the pulses in 
a pulse pair can be adjusted by a delay line. For the best 
signal to noise contrast in the microscopy, we set r close 
to zero. The electric field after another beam splitter 
that combines the two pulses collinearly is given by 

Eiot{t) = ^Mt){cOs{{LO + + COs((w + ( 1 ) 
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An n-photon absorption is proportional to (F^tot(^))^’^ 
(the complete expression for n-photon absorption is given 
in the Appendix). The intensity of the resulting flu¬ 
orescence gets modulated at the frequencies given by 
/mod = xpi ± yp 2 ^ where x and y are non-negative in¬ 
tegers less or equal to n, i.e. x G {0,1, 2,..., n} and 
y G {0,1, 2,..., n}. We set pi = 55 MHz and p 2 = 55.1 
MHz in our experiments. The avalanche photodiode 
(APD) used in the setup has a bandwidth of 2 MHz. Con¬ 
sequently only the fluorescence modulated at the multi¬ 
ples of the relative frequency (A0 = p 2 — Pi) can be 
resolved in our measurements. The amplitude as well as 
the phase of the signal can be obtained using a signal de¬ 
modulator such as a lock-in amplifier. We use generalized 
lock-in amplifiers for this purpose [TTHT3]. 



FIG. 1. Schematic of the experimental setup. In the fig¬ 
ure OS: Oscillator, BS: Beam splitter, M: Mirror, AOM: 
Acousto-optic modulator, DL: Delay line, PD: Photodiode, 
MICR: Microscope, FIL: Band pass filter (550± 25 nm), 
APD: Avalanche ohotodiode, SIG: Signal, REF: Reference 
and SDM: Signal demodulator. 

The modulated intensity of the two photon fluores¬ 
cence recorded in the measurements is given by 

Jdet(i) oc (4cos((A^)i) + cos(2(A^)i)). (2) 

The fluorescence is separated from the excitation beams 
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by using a band pass filter with pass band at 550± 25 nm. 
The center wavelength of the excitation beams is at 780 
nm and the full width at half maximum of the spectrum 
is about 130 nm. Figure [^a) shows the Fourier trans¬ 
form of the fluorescence of CdSe quantum dots recorded 
for 10 ms. The transform shows two peaks at 0.1 and 
0.2 MHz, which correspond to the contributions at A0 
and 2A0 frequencies, respectively. The two contributions 
arise from the two different field matter interaction path¬ 
ways. Out of the four field interactions when two of the 
fields come from each of the beams, the fluorescence gets 
modulated at 2A0, when three of the fields come from 
one of the beams and the fourth from the other beam the 
fluorescence gets modulated at A0 (see the Appendix for 
detailed discussion). Equation predicts that the con¬ 
tribution at A0 is four times the contribution at 2A0. 
The amplitude of the fluorescence at the two frequencies 
in the measurements (Figure |^a)) deviate slightly from 
this relation. Equation is valid when the intensity of 
both the beams at the focus spot are equal. However, 
this is not always the case. In general the signal at A(/) 
is more than or equal to four times the signal at 2A0 
(see the Appendix for further discussion). The signals at 
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FIG. 2. (a) Fourier transform of the two photon fluorescence 
from CdSe quantum dots recorded for 10 ms. The phases of 
the two beams used to excite the fluorescence are modulated 
at 55 and 55.1 MHz, respectively. The resulting fluorescence 
is modulated at 0.1 and 0.2 MHz. (b) noise of the detector at 
different frequencies, (c) raw image of a cluster of CdSe QDs 
taken by recording the normal two photon fluorescence image 
and (d) image of the same cluster taken by using the phase 
modulation of the beams. 

both the frequencies, A0 and 2A0, can be used in mi¬ 
croscopy. Figure [^(c) shows the raw data of the normal 
two photon fluorescence image of a cluster of CdSe quan¬ 
tum dots on a glass slide. The image is taken by raster 
scanning the focus spot over the sample. The dwell time 


per scan point is 10 /is. We use a reflective objective 
(Adjustable ReflX Objective 36X/0.5NA IR, Finite, Ed¬ 
mund optics) mounted on an inverted microscope (Nikon 
Ti-S) to focus the laser beams onto the sample. The 
focal length of the objective is about 5 mm. The fluo¬ 
rescence is collected in the epi-direction. A multiphoton 
short-pass dichroic beamsplitter with the edge at 670 nm 
(EE670-SDi01-25x36, Semrock) separates the excitation 
beam from the collected fluorescence. The band-pass fil¬ 
ter at 550 nm (OD 4) before the detector further reduces 
the chance of excitation photons reaching the detector. 
Eigure (d) shows the raw image of the same cluster 
using phase modulation. The image is recorded at the 
modulation frequency of 0.1 MHz. Comparing the two 
images, we see that the image taken with phase mod¬ 
ulation has reduced noise.The reduction in the noise is 
mainly due to the fact that at low frequencies the sig¬ 
nal is heavily contaminated by the so called “pink” noise 
(1// noise) that is prevalent in all the electronic systems. 
Eigure]^ (b) shows the electronic noise of the signal de¬ 
tection system that include the APD, the cables and the 
digitizer at different frequencies. 

Apart from the improvement in the signal to noise 
ratio, a major advantage of using phase modulation in 
MM is in multiplexing the raster scanning by collect¬ 
ing fluorescence from multiple excitation spots in the 
sample. Various techniques have been developed for 
multi-focal multiphoton microscopy (MMM) over the last 
decade El ig 0 d- The different techniques have their 
strengths, such as high speed imaging m, and weak¬ 
nesses, such as inhibited depth scanning [8]. The most 
recent advances include temporal multiplexing of the sig¬ 
nal [8] to reduce the effect of interferences between the 
spatially overlapping foci, and remote focusing [151 [16] fo 
allow rapid axial scanning. When applying phase mod¬ 
ulation in MMM, we use more than one pair of beams 
to illuminate different spots on the sample. Each pair 
of the beams have different relative modulation frequen¬ 
cies, A(/)i, A(/) 2 , .... The fluorescence from all the spots 
are directed to a single element detector such as a photo¬ 
multiplier tube or an APD. The use of single element de¬ 
tector eliminates the limitations on the axial resolution 
which is inherent in imaging techniques that use array de¬ 
tectors Bm- The signals from the different spots can be 
isolated by demodulating them at their respective mod¬ 
ulation frequencies. This method of multiplexing signal 
in the frequency domain has advantages over the meth¬ 
ods that use multiplexing in the time domain. Usually 
for the temporal multiplexing, the time between the dif¬ 
ferent beams need to be delayed by the response time 
of the detector, which in most of the cases is a couple of 
nano-seconds |8] . This requires a rather elaborate optical 
setup or a specialized light source. Moreover, only few 
beams can be multiplexed within a repetition period of 
about 10-12 ns of an oscillator. 

Eigure shows two images of a cluster of CdSe quan¬ 
tum dots. The two images are produced by raster scan¬ 
ning two pairs of beams. The relative phase modulations 
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are set to A(/)i =0.1 and A (/)2 = 0.17 MHz, respec¬ 
tively. Only the images acquired at 2A0i = 0.2 (top) 
and 2A02 = 0-34 MHz (bottom) are shown. The two 
pairs of the beams are focused at different axial position 
(depth) separated by 4 /im on the sample. The distance 
between the focal points is controlled by remote focusing 
(i.e. by changing the divergence of one of the pairs of the 
beams with respect to the other). 



FIG. 3. Multi-point raster scanned images of a cluster of 
CdSe quantum dots on a glass slide. Image (b) is acquired at 
a depth of 4 /im below image (a). The relative phase modu¬ 
lations are set to 0.1 and 0.17 MHz, respectively. The images 
themselves are acquired at twice the relative modulation fre¬ 
quencies. 

There is a technical difference between the normal and 
the phase modulated MMM in the way information is 
acquired from each raster scanned point. In a normal 
MMM one avoids any spatial and temporal overlap of 
the different beams m- In the phase modulated MM 
one gets the information from the overlapping regions. 
This difference could be exploited to increase the num¬ 
ber of points from which the signal is sampled. Instead 
of focusing two pairs of beams (four beams altogether) 
to two separate foci one can also use a 2D beam pat¬ 
tern as shown in Figure [^a) where the signal from four 
overlapping regions can be recorded at once. In general 
a square grid of m x n beams can have 2mn — m — n 
overlapping regions. The number of over-lapping regions 
can be increased further by using a 3D grid. Figure [^b) 
shows a 2 X 2 X 2 grid of beams with 12 overlapping re¬ 
gions. In general a cubic grid of m x n x o beams can have 
3mno — mn — no — om overlapping regions that can be 
used in multi-point scanning in phase modulated MMM. 

To conclude, we have shown how the phase modulation 
of the laser beams of ultra-short pulses can be used to 
modulate the intensity of the non-linear interaction with 
the sample. We have also shown that phase modulation 
can be used in MM to improve the signal to noise ratio 
and to multiplex the signals in the frequency domain for 
MMM. Though we have demonstrated the principle using 
two photon fluorescence, the phase modulation can be 
applied to any non-linear signal, such as second harmonic 
generation, to achieve similar advantages. 
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FIG. 4. A two dimensional grid (a) and a three dimensional 
grid (b) of beam patterns with over-lapping regions that can 
be used in multi-point raster scanning, (a) has four overlap¬ 
ping regions and (b) has 12 overlapping regions. The number 
of overlapping regions increases with the dimensionality of the 
beam pattern. 
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Appendix: Expression for multi-photon ab¬ 

sorption using two phase modulated beams. Let 

(cj + 0i) t = a and {uj (j) 2 )t = b. The expression for 
n-photon absorption is given by: 

^to"t it) = it) (^2 cos (^) 

= (t) cos^^ (x) cos^^ (y) 

with X = (a + 6) /2,// = (a — b) /2. Using the trigono¬ 
metric identities for the n-th power of cosine functions, 
we get: 

/ n — 1 

sS it) oc ( + 2 ^ Cf cos {{2n - 2k) x) 

k =0 
n — 1 

+ 2 cos ((2n - 21) y) 

1=0 

n — 1 

= (Cn”)^ + 2C'^’" cl"' cos ((2n - 2k) x) + 

/c=0 

n — 1 

2C'^’" y] cl"' cos ((2n - 21) y) 

k =0 

. n—1n—1 

+ ^ T cos ((2n - 2k) x) cos ((2n - 21) y). 

k =0 1=0 

Only the fluorescence arising from the third term (term 
that has (2n — 2l)y as the argument of the only cosine) 
in the expression given above contains modulations that 
can be resolved by a slow detector. Reverting back to 
the original notations, we get 

n—1 

Idet it) OC y] Cf” COS ((n - 1) {(p2 -(pl)t)- 
/=0 

The expression simplifies to 

Idet it) CC cos {2 {4>2 - 4>l)t) +4 cos {{4>2 - 4>l)t) 
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for the fluorescence due to two photon absorption. 

When the intensities of the two beams are not the same 
we can follow the similar derivation with different en¬ 
velopes of the fields, Ai (t) and A 2 (t). The expression 
for the resulting fluorescence intensity is given by 

Idet {t) oc Aj {t) Al {t) cos (2 (02 -4>l)t) + 

2Ai (t) Al (t) cos ((02 - 0i) i) + 

2Al (t) A 2 (t) cos ((02 - 0 i) i), 

which shows that when three of the fields interacting with 
sample are from one of the beams and the fourth field is 
from the other beam, the two photon fluorescence gets 
modulated at 02 — 01- Similarly, when two of the fields 
from each of the beams participate in the non-linear ab¬ 


sorption, the fluorescence gets modulated at 2(02 — 0i). 
If we take the general identity 

{A^{t)-A2{t)f>0 

then we get 

^l(^) ^2(^) — ' 

As Ai{t) and A 2 {t) are non-negative real numbers, mul¬ 
tiplying both the sides by 2 Ai{t)A 2 {t) gives 

2Al{t)A2it) + 2Ai{t)Al{t) > 4Al{t)Al{t). 

So in general, the fluorescence modulated at 02 — 0i is 
more or equal to four times the fluorescence modulated 
at 2(02 - 0 i). 
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